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We consider the problem of electron transport in segregated conductor-insulator composites in 
which the conducting particles are connected to all others via tunneling conductances thus forming 
a global tunneling connected resistor network. Segregation is induced by the presence of large 
insulating particles which forbid the much smaller conducting fillers to occupy uniformly the three 
dimensional volume of the composite. By considering both colloidal-like and granular-like dispersions 
of the conducting phase, modeled respectively by dispersions in the continuum and in the lattice, we 
evaluate by Monte Carlo simulations the effect of segregation on the composite conductivity a, and 
show that an effective medium theory applied to the tunneling network reproduces accurately the 
Monte Carlo results. The theory clarifies that the main effect of segregation in the continuum is that 
of reducing the mean inter-particle distances, leading to a strong enhancement of the conductivity. 
In the lattice segregation case the conductivity enhancement is instead given by the lowering of 
the percolation thresholds for first and beyond-first nearest neighbors. Our results generalize to 
segregated composites the tunneling-based description of both the percolation and hopping regimes 
introduced previously for homogeneous disordered systems. 

PACS numbers: 64.60. ah, 73.40.Gk, 72.80.Tm, 72.20.Fr 



I. INTRODUCTION 

Transport properties of conductor-insulator compos- 
ites are strongly influenced by the microstructural char- 
acteristics of the composite itself. In particular, the dc 
conductivity a depends on the volume fraction </> of the 
conducting constituents, on their size [l| and shape 
as well as on their spatial arrangement in the composite 
[1, 0] • All those aspects influence a through their effects 
on the electrical connectivity of the conductive phase, 
and can thus be exploited to meet speciflc criteria for 
the transport properties in composites. 

In compacted mixtures of micrometric conducting and 
insulating powders the electrical connectivity is es- 
tablished by direct contact connections between the con- 
ducting particles In this situation, a displays an 
insulator-conductor transition when the concentration (p 
of the conducting phase is such that a macroscopic cluster 
of connected particles spans the entire sample, allowing 
the charge carriers to flow between the electrodes. Perco- 
lation theory 0, describes such transition by mapping 
the inter-particle electrical connections to a random re- 
sistor network, where the elemental conductances g are 
either when there is no contact oi g ^ Q when two 
conducting particles touch each other. In this way, per- 
colation theory predicts a power-law behavior of the form 
(T ~ ((/) — (pcY for (j) > (f)c^ where (pc is the critical volume 
fraction beyond which a spanning cluster is formed and t 
is a universal transport exponent taking the value i ~ 2 
for all three-dimensional systems For this kind of 

composites, a is thus mainly controlled by the value of 
<f)cj which depends upon the shape and the dispersion of 



'Electronic address: 
t Electronic address: 



biagio.mgro@epfl.ch| 



claudio.grimaIdi@epfl.ch| 



the conducting particles. 

In nanocomposites made of colloidal dispersions of 
nanometric conducting particles in an insulating contin- 
uous matrix the predominant mechanism of transport is 
not by direct contact, but rather through indirect elec- 
trical connections between particles established by quan- 
tum tunneling [llj . For temperatures sufficiently high to 
neglect particle charging effects, energetic disorder, and 
Coulomb interaction, the inter-particle conductance de- 
cays exponentially for distances between particles larger 
than a characteristic tunneling length ^, which measures 
the electron wave function decay within the insulating 
phase. Although ^ depends on specific material proper- 
ties, its value is nevertheless limited to a few nanometers 
or less, which is a relevant length scale for composites 
with nanometric conducting particles, whose typical sizes 
range from tens up to hundreds of nanometers. Hence, 
besides the effect of shape and dispersion of the conduct- 
ing constituents, the composite conductivity a is in this 
case affected also by the mean distance between the par- 
ticles, which essentially depends on (j). 

It is clear that tunneling conduction, albeit decaying 
fast, does not imply a sharp cutoff of the connectivity 
between particles and the introduction of a contact-like 
connectivity criterion, as done for powder mixtures, is 
not suitable [l2j. Conversely, inter-particle conductance 
by contact can be seen as a limiting case of tunneling 
when the particle sizes D are much larger than the tun- 
neling decay length ^, as it is the case for mixtures of 
micrometric conducting and insulating powders [^-@. It 
turns out that it is indeed possible to describe both con- 
tact and distance dependent connectivity mechanisms 
within a single formalism, in which the conducting parti- 
cles are all electrically connected to each other by tunnel- 
ing [ll]. By using this global tunneling network (GTN) 
approach, percolation properties of compacted powders 
or of other granular materials with micrometric conduct- 



ing grain sizes can then be recovered by the D/^ ^ 1 
limit of the theory, while hopping transport of colloidal 
nanocomposites is obtained by much smaller D/S^ values. 

As noticed in Ref. [ll], D/^ is not, however, the only 
factor discriminating between percolation and hopping 
regimes. Indeed, for D/^ sufficiently large, the com- 
posite microstructure plays the most relevant role and, 
through the arrangement of the conducting matter in the 
composite, promotes one regime or the other. For exam- 
ple, homogeneous dispersions of impenetrable conducting 
spheres in the (insulating) continuum are always charac- 
terized by an hopping (or, cquivalcntly, tunneling) type of 
transport. In this case indeed the conductivity decreases 
fast but continuously as 4> is reduced because, in average, 
the inter-sphere distances increase. Instead, percolation- 
like behavior of transport arises in close-packed mixtures 
of conducting and insulating spheres because contact or 
near-contact clusters of conducting spheres span the en- 
tire system for all volume fractions larger than (pc- For 
large D/S^, multiple percolation thresholds due to clusters 
of further next-nearest neighbors can arise in fractionally 
occupied periodic lattices. 

In this article we extend the study of the percola- 
tion and hopping regimes to the case of (locally) non- 
homogeneous dispersions of conducting particles by con- 
sidering segregation of the conducting fillers due to large 
(compared to the conducting particles size) insulating in- 
clusions [l^ - [l7| . This type of microstructure is rather 
common in real composites and is at the origin of the 
large conductivity values measured also for very low con- 
tents of the conducting phase. We shall show how this 
peculiar behavior arises in both hopping and percolat- 
ing regimes by solving numerically the tunneling network 
equations for continuum and lattice segregated particle 
distributions. Furthermore, by using a generalization 
[H, [13] of the classical effective medium approximation 
(EMA) O, m, [131, we explicitly relate the microstruc- 
ture properties of the composites with the transport be- 
havior, and provide an approximate but accurate analytic 
treatment of the conductivity problem, thus extending 
our previous results for the continuum jlTf and general- 
izing the formulation to the segregated case. 

In Sec. |TT] we consider the continuum regime by de- 
scribing the model for the segregated composites and by 
introducing the EMA formulation for the calculation of 
the overall conductance. In Subsec. Ill Al we present the 
results of both EMA calculation and MC simulations for 
some segregated systems and in Subsec. IlIBI wc provide 
an explicit approximate analytical formula for the com- 
posite conductivity based on EMA. Section. |lll]is devoted 
to the calculation of both EMA and MC for lattice seg- 
regated dispersions of conducting fillers. Finally, Sec. IIVI 
is devoted to the conclusions. 
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FIG. 1: (Color online) Schematic two-dimensional representa- 
tion of the segregated dispersions in tlie continuum (a) and in 
the lattice (b) . The small black circles denote the conducting 
spheres of diameter D\ while the larger penetrable circles are 
the insulating spherical inclusions of diameter D2. In (b) the 
centers of the conducting spheres occupy a fraction of the site 
of a periodic lattice. 



II. SEGREGATION IN THE CONTINUUM 

Segregation is very common in composite materials 
such as Ru02-based cermets [2l| or polymer-based com- 
posites This effect is particularly evident when the 
mean size of the insulating grains (glass-particles or poly- 
meric inclusions as large as few micrometers) is much 
larger than that of the conducting fillers, whose typical 
size ranges from tens to hundreds of nanometers. The 
presence of such large insulating inclusions reduces the 
volume available for the conducting fillers, which are con- 
fined in the remaining space, leading thus to a locally 
non-homogeneous distribution of the conducting phase 
in the composite. 

In the following we represent the conducting phase of a 
conductor-insulator composite by generating dispersions 
of Ni hard-sphere particles of diameter Di and volume 
fraction 0i — pivi, where vi = tt 13^/6 is the volume of a 
single sphere, pi — Ni/V is the number density, and V is 
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the total volume. The non-segregated (or homogeneous 
case) is obtained through a random sequential addition 
(RSA) of hard spheres into a cubic box of side L, fol- 
lowed by an equilibration Monte Carlo (MC) process as 
described in Ref. [T^. When the RSA limit ^^^^'^ - 0.38 
is reached [2^, the equilibrium configuration is obtained 
by placing initially the hard spheres into a cubic lattice 
and then by relaxing the system through MC runs. 

The segregated regime is schematically shown in 
Fig.[lja) and is composed by a mixed system of mutually 
impenetrable conducting and insulating spherical parti- 
cles. This is obtained by considering a random dispersion 
of N2 fully penetrable spheres of diameter D2 represent- 
ing the insulating inclusions [IB, ■ We fix the diameter 
ratio of the two particle species at D2/D1 = 8, having no- 
ticed that this condition is sufficient to characterize 
segregated systems in the D2 ^ -Di regime. In addition, 
in order to minimize size effects, L is chosen as to be at 
least one order of magnitude larger than D2- Since the in- 
sulating inclusions are placed randomly and are penetra- 
ble, the occupied volume fraction is </)2 = 1 — exp(— i;2P2) 
[H, where V2 = and p2 = N2/L^. 

After having placed the insulating spheres, the remain- 
ing available space is filled with given densities pi of 
conducting hard spheres by means of the same place- 
ment and equilibration procedures as for the homoge- 
neous case. Given the mutual impenetrability of the two 
kinds of spheres, the available volume fraction for ar- 
ranging the centers of the conducting fillers is given by 
V* = exp(-T4xP2), where = 7i'(-Di -I- £'2)^/6 is the 
excluded volume of an insulating sphere with respect to 
a conducting one. By exploiting the previous definition 
of 62, V* thus reduces to 1161: 



(1 



\(l+Dt/D2f 



(1) 



which defines an effective filler volume fraction (pcs = 
(j)ilv*. 

As stated in the introduction, we consider all con- 
ducting particles as electrically connected to all others 
through tunneling conductances which, for two generic 
impenetrable spheres i and j of diameter Di placed at 
positions and r^, is given by: 



9ij = go exp 



(2) 



where go is a constant "contact" conductance which will 
be set equal to the unity in the following, ^ is the charac- 
teristic tunneling decay length, and r.y = jr^ — r^l is the 
distance between two conducting sphere centers. 

Contrary to classical resistor networks with few con- 
nected nearest neighbors [l^ , the ensemble of all tunnel- 
ing conductances forms a fully connected weighted net- 
work of A^i nodes, each having A^i — 1 neighbors. The 
overall conductivity depends on Di and ^, as well as 
on the volume fraction (jji of the particular conducting 
sphere distribution. Note that the aforementioned map- 
ping between the sphere system and the resistor network 



holds for both the homogeneous and segregated disper- 
sions of the conducting fillers since the information about 
the specific distribution is implicit in the weighted links 
[i.e., the conductances gij of Eq. ^] over all sphere cen- 
ters. 



A. EMA and MC 

All these dependencies, and in particular the relation- 
ship between the spatial particle arrangements and the 
global transport properties, can be made explicit by em- 
ploying the EMA formulation developed in Refs. [T3l.lT8|. 
which is a generalization to complete tunnelin g re sistor 
networks of the classical EMA approach [lol. [l9l.l20| . The 
original tunneling network is thus replaced by an effective 
one where all bond conductances are equal to g, whose 
value is found by requiring that the effective network has 
the same average resistance as the original. By consider- 
ing only two-site clusters [l^, the following equation for 
g is found: 



9ij - 9 



+ [{N,-l)/2-l]d 



0. 



(3) 



where (. . .) indicates a configurational average. By in- 
troducing the radial distribution function (RDF) (72 (j') 
of the Ni conducting hard spheres defined as [25| : 



P92{r) 



equation ([3]) can be recast as follows 



dr 



Anr^pi g2{i') 



5*cxp[2(r-i?i)/e] + l 



(4) 



(5) 



where g* ~ Nig/2 is the conductance between any two 
nodes of the effective network. According to Eq. ([5]) , the 
whole information about the spatial distribution of the 
conducting particles is contained in g2 {r) , which therefore 
directly governs the behavior of the EMA conductance 
g* . To illustrate how the segregation affects g* through 
its effects on 52 (t), we plot in Fig. [2] the numerically cal- 
culated RDF for = 0.114 and for three different values 
of the volume fraction (j)2 of the insulating spherical in- 
clusions. The (/)2 = case corresponds to a homogeneous 
dispersion of hard spheres in the continuum, and the re- 
sulting RDF (circles) is basically featureless for r > Di 
(apart for a slight increase near contact) because of the 
rather low value of 0i used. For 02 > two features 
emerge. First, the RDF develops an oscillating behavior 
with a clear peak at r ~ 2Di and a second one, visible 
for 02 = 0.6, at r ~ 3Di. Second, the RDF is enhanced 
with respect to the homogeneous case for all values of 
r lower than r ~ D2 = 8D1. These characteristics re- 
semble in part those observed in the RDF of the smaller 
particles in hard-core mixtures [26| and are indications 
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FIG. 2: (Color online) Numerical RDF as a function of the 
distance in the continuum regime for three different volume 
fractions of the insulating phase (j>2 ~ 0, 0.4, and 0.6, for 
D2/D1 — 8, and for a fixed filler volume fraction (jii — 0.114 



of the enhanced probabiUty of having particles separated 
by a distance r < D2 when 02 > 0, which means that 
the conducting fillers are in average closer to each other 
when segregation increases. Due to the exponential de- 
pendence on the particle separation of the tunneling 
conductance, Eq. ([2]), the EMA conductance at fixed (pi 
is thus expected to be increased by segregation because, 
for 02 > 0, particles have lower average value of r.y. 

This is confirmed by the results plotted in Fig. ^a), 
where the EMA conductance g* , obtained from solving 
Eq. ([5]) with g2ir) evaluated from MC calculations, is 
reported as a function of 0i and for 02 = 0, 0.4, and 0.6. 
The role of segregation in enhancing g* is made even 
more evident when Di/^ increases, as it can be inferred 
by comparing the results in the upper panel of Fig. El^a), 
obtained for Di/^ = 15, with those with Di/^ = 50 in 
the lower panel. 

The result that the EMA formulation provides a trans- 
parent relation between the microstructure of segregated 
continuum composites, contained in g2{f'), and the trans- 
port behavior is made even more firm by the fact that 
the EMA conductance is in excellent accord with our 
fully numerical calculations of a [13]. These are shown 
in Fig.[3ljb), and have been obtained by the same numeri- 
cal procedures described in Ref. (l7| consisting in solving 
numerically the Kirchoff equations of the tunneling re- 
sistor network with conductances given by Eq. O For 
the segregated cases, each symbol in Fig. [31[b) is the out- 
come of Njf = 200 realizations for A^i conducting parti- 
cles ranging from a few hundreds for low (pi to hundreds 
of thousands for the higher values, by keeping L fixed. 
Instead for the homogeneous case (02 ~ 0) ~ 500 
has been considered, with A^i fixed at 2± 1000. 



B. EMA analytical formula 

The good agreement between the EMA results and the 
fully numerical cr in both the homogeneous and segre- 
gated regimes suggests that further insights can be gained 
directly from the EMA equation ([5]). Hence we proceed 
here with some further approximations in the attempt to 
find an analytical expression for the EMA conductance. 
We start by noticing that the integral in Eq. ([5]) can be 
rewritten as 



drinr^ pi g2{'>')W{r) 



where 



W{r) 



exp 



(6) 



(7) 



and r* is defined by the following relation: 



exp 



(r* - Di) 



(8) 



We note that for large values of Di/^, which is the regime 
of practical interest for our purposes [l3| (see also SecH]), 
W{r) is well approximated by 9{r*—r), where 6{x) is 1 for 
X > and otherwise. Thus, by adopting the definition 
of the cumulative coordination number 



Zir) = / dr'4^r'Vi52(r'), 



(9) 



which gives the number of spheres whose centers are 
within a distance r from the center of a given sphere, 
we obtain the following approximation of Eq. ([5|) : 



Z{r*) = 2. 



(10) 



Equation (|5]), where r* is such that Eq. pUj) is satisfied, 
is plotted in Fig.[3l^a) (star symbols) and is in close agree- 
ment with the numerical solution of Eq. ([5]). We proceed 
further by noticing that the quantity 6* = r* — Di is 
the EMA equivalent of the critical distance Sc found by 
the critical path approximation (CPA) for the tunnel- 
ing conductivity [ll,!!!]. The CPA 4, plotted in Fig. g] 
(symbols) as a function of 0i for 02 = 0, 0.4 and 0.6, 
is the shortest among the interparticle distances such 
that the subnetwork formed by those particles having 
rij — Di < Sc forms a percolating cluster. Equivalently, 
dc is such that Z{Sc + Di) = Z^ is satisfied, where Z^ is 
the critical coordination number which, for the homoge- 
neous case 02 = 0, ranges between Zc ~ 2.7 for 0i 
and Zc — 1.5 for 0i ~ 0.5 [13 • It turns out therefore 
that the right-hand side of Eq. ([TU]) falls well within the 
range of possible Zc values, which is the ultimate reason 
of the good accord between the EMA approximation of 
the conductance and the fully numerical a already no- 
ticed in Ref. [l^ for the homogeneous case. 
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FIG. 3: (Color online) (a) Calculated EMA conductance g* (open symbols) as a function of the volume fraction 01 of the 
conducting spheres of diameter Di for -Di/^ = 15 and 50, and for different values of the volume fraction 02 = 0, 0.4, and 0.6 
of the insulating spheres with diameter D2 and diameter ratio D2/D1 = 8. The crosses and the dashed lines refer respectively 
to Eq. (IS}, where r* is such that Z{r*) = 2 is satisfied, and to Eq. p4|l . (b): conductivity o for the same parameters as in (a) 
obtained from the numerical solution of the tunneling resistor equations. 



By noticing from Fig. |4] that 5* <^ Di for large (f>i 
values and that g2ir) = for r < Di, in order to cor- 
rectly capture the high density regime, we approximate 
the RDF in Eq. ^ with its contact value 92(^1)- In this 
way, Eq. ([T0| can be rewritten as 



92{Di) drWpi%-^i) =2, (11) 

Jo 



which can be solved for 5* leading to: 



-,1/3 



-A 



(12) 



By using the Carnahan-Starling formula (?2£Di) = (1 — 
(/)i/2)/(l - 0>i)3 for the RDF at contact [3l|, which is 
a well-known approximation used in the theory of sim- 
ple liquids [1^, Eq. ([T^ turns out to be a rather good 
approximation for dc in the whole range of densities for 
(j)2 = (solid line in Fig. 2]), and thus we can use the 
same approximate scaling relation that we formulated in 
Ref. [13 for the critical distance Sc of segregated systems. 



Hence, if (5*(0i, v*) is the EMA critical distance for a seg- 
regated system parametrized by the available volume v* 
of Eq. (P, then 



') ^ ^;*-l/3^*(0,ff), 



(13) 



where (f>cs = 0i/"* is the effective volume fraction for 
the conducting fillers introduced in Sec.|lTl Equation 
then states that the EMA critical distance in the segre- 
gated regime can be directly obtained from that of the 
homogeneous case calculated at (j)cs. As shown in Fig. 21 
Eq. (|13p compares well with MC calculations of the crit- 
ical distance 6c, so that by using Eq. with Eqs. 
and (|13p and 92(^1) as given by the Carnahan-Starling 
expression we obtain the following approximated analyt- 
ical formula of the EMA conductance: 



g = exp • 



1 + 0cff + 



20off (2 - 0c 



- 1 



(14) 

The above expression is plotted in Fig. [DJa) by dashed 
lines and turns out to be in very good agreement with 
the full solution of the EMA integral of Eq. ([5]). Hence, 
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FIG. 4: (Color online) Critical distance Sc dependence on the 
volume fraction cpi of the conducting spheres for (f)2 = 0, 0.4, 
and 0.6 and for D2/D1 = 8 extracted from MC calculations 
[l7| . Solid lines: our 5* approximation from Eq. (|12p (02 ~ 0) 
and from Eq. ^ (02 / 0). 



despite of its simplicity, Eq. ([14]) effectively captures 
the conductance behavior of both the homogeneous and 
segregated cases in the whole r ang e of 0i values, thus 
generalizing the results of Refs. |l2l. [isl. [l7j. 

III. SEGREGATION IN THE LATTICE 

Segregation can be induced in powder mixtures of con- 
ducting and insulating particles when the mean size of 
the insulating grains is much larger than that of the con- 
ducting ones. Typical examples of segregated powder 
mixtures and of the corresponding conductivity measure- 
ments can be found in Refs. [H, [J. Compared to powder 
mixtures of conducting and insulating grains with com- 
parable mean sizes @,13|: the conductivity of segregated 
powders drops by several orders of magnitude at much 
lower values of the volume fraction. 

In order to describe the effect of segregation in pow- 
der mixtures, we consider a model composite where the 
conducting particles occupy only a fraction p of the total 
M sites of a simple cubic lattice. For non-segregated sys- 
tems, we consider equal sized conducting and insulating 
particles of spherical shape with diameter equal to the 
lattice spacing, as in the Scher and Zallen model [s^ . 
For a random distribution of the conducting spheres on 
the cubic lattice the corresponding RDF reduces to [l3| : 

/'•92M = £ E ^Sir-R,), (15) 

k=l,2.... * 

where A/fc is the number of the fcth nearest neighbors 
being at distance Rk from a reference particle set at the 
origin. From Eq. (|15p . and by using the tunneling inter- 
particle conductance of Eq. the EMA equation ([S]) 



becomes 

^.5....5*e^P[2(^fe-^i)/e] + l ^^^^ 

whose solution is plotted in Fig. [SJa) (solid black lines) 
as a function of the volume fraction (jji of the conducting 
phase (01 = pir/G) for two different values of -Di/C- As 
discussed in more details in Ref. [l3|, the decrease of g* 
as 01 is characterized by sharp drops at 0i = pfe7r/6 
where pk = 2/(^^,^-|^ A/^') is the percolation threshold 
for the fcth nearest neighbors. This behavior is due to 
the discrete nature of the lattice RDF, Eq. (jTS]), and 
is confirmed by the full MC results shown in Fig. ^h), 
where the conductivity is obtained by following the same 
procedure of the continuum regime for the number of 
particles held fixed at A^i ^ 1000. 

Since in real segregated powder mixtures the conduct- 
ing particle sizes are in the micro- metric range 
the corresponding large values of -Di/^ make the perco- 
lation threshold for particles at contact pi = 2/A/i (= 
1/3 for a cubic lattice) the one of practical interest. In 
this regime, and for 0i > 0j , the EMA conductance and 
the Monte Carlo conductivity follow the power-law be- 
haviors g* oc (01 — 0}) and a (x (0i — (f>\Y, with respec- 
tively 0j = 7r/18 ~ 0.174 and 0} ~ 0.163, where t ~ 2 
is the transport exponent for three-dimensional systems 

In analogy with the continuum segregation of Sec. [ill 
we consider the segregation in the lattice as being due 
to 7V2 penetrable and insulating spheres of diameter D2 
placed at random in the cubic volume. As schematically 
shown in Fig.[TJb), the conducting particles will thus oc- 
cupy randomly only those lattice sites lying outside the 
insulating spheres and not leading to overlaps between 
the two species of particles. As for the homogeneous 
random lattice case, the RDF of the segregated lattice 
is given by a series of delta-peaks centered at Rk, as in 
Eq. ([15]), but with the number A/fc of the fcth nearest 
neighbors being dependent of D2 /Di and of the volume 
fraction 02 of the insulating spheres. In Fig. [5] we show 
the Rk dependence of A/fc (02) in units of the number 
of the fcth nearest neighbors for the homogeneous lat- 
tice A4(0) for D2/D1 8 and for 02 = 0.4 and 0.6. 
For Rk » Di the ratio A/'fc(02)/A/fc(O) approaches unity 
which indicates that at large distances segregation plays a 
minor role, a result similar to the one found for the con- 
tinuum segregation case (see Fig. [2]). On the contrary, 
for Rk values close to contact, segregation enhances the 
number of fcth nearest neighbors for fixed 0i because, 
again in analogy with the continuum case, the reduced 
available volume enhances the probability of finding con- 
ducting particles at distances lower than about D2. By 
presuming that the enhanced local probability can be 
approximated by Pi =pi/v*, then A/fe(02)/A/fc(O) should 
scale as 1/v* for Rk sufficiently close to contact (and 
D2/D1 sufficiently large), which is a fair approximation 
for 02 = 0.4 but a less satisfactory one for 02 = 0.6 
(horizontal dashed lines in Fig. [H]). 
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FIG. 5: (Color online) (a) Calculated EMA conductance g* for the segregated lattice as a function of the volume fraction (t>i 
of the conducting spheres for -Di/^ = 15 and 50, and for different values of the volume fraction 02 of the insulating spheres 
with diameter D2. (b): conductivity a for the same parameters as in (a) obtained from the numerical solution of the tunneling 
resistor equations. 



From pkXh) - 2/Et^iA4'(<^2)] and Afk{(f>2) > 
7\4(0) it follows immediately that the percolation thresh- 
olds Pki(t>2) tor the fcth nearest neighbors in the segre- 
gated lattice are lower than those for the homogeneously 
random lattice Pfc(O). In particular, the percolation 
threshold of particles at contact is reduced by the fac- 
tor A/'i(O)/7Vi(02) which from Fig. [His 0.54 and 0.33 for 
(j)2 = 0.4 and (/)2 = 0.6, respectively. Note that approx- 
imating p\{(l)2)/p\iO) with V* leads to 0.48 for (j)2 — 0.4 
and to 0.27 for 02 = 0.6. 

The systematic lowering of Pk{4'2) as (j)2 is enhanced is 
reflected in the (/)i-dependence of the EMA conductance 
g* in Fig. obtained by solving numerically Eq. 
with our calculated values of N'ki(f>2)- The overall ef- 
fect of segregation is thus the enhancement of 5* with 
respect to the homogeneous lattice case at 02 = in- 
duced by the downshift of all Pk{4>2) values. This behav- 
ior is confirmed by the full MC results of Fig.[Sl[b), which 
have been obtained by solving the tunneling resistor net- 
work for Ni ranging from a few hundreds for 0i ~ 10""^ 
to Nx - 230000 for (px/v* - 0.5 with L held fixed at 
L = 10-D2- Given the above results and discussion, the 
conductivity for segregated micrometric (i.e., Di/^ ^ 1) 



powders just above the percolation threshold is expected 
thus to follow approximately a oc (0i — (f)\v*Y where 
0} is the critical volume fraction for conducting parti- 
cles at contact in the absence of segregation. Due to the 
quasi-invariance of 0}, according to which 0} « 0.17 in- 
dependently of the (three dimensional) lattice topology 
j32| . this result is expected to apply also to non-cubic 
segregated lattices. 



IV. CONCLUSIONS 

In this paper we have considered the dc electrical trans- 
port problem in two-phase segregated amorphous solids, 
where large insulating inclusions prevent the smaller con- 
ducting particles to be dispersed homogeneously in the 
three dimensional volume. By taking into account ex- 
plicitly the tunneling mechanism of electron transfer be- 
tween conducting particles, we have studied the effect 
of segregation for both continuum and lattice models of 
composites. For continuum segregated composite mate- 
rials, we have shown by theory and Monte Carlo simula- 
tions that segregation basically reduces the inter-particle 
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FIG. 6: (Color online) Enhancement ratio of the number TVfe 
of fcth neighbors for 02 = 0.4 and 0.6 (for D2/D1 = 8) as 
compared to the non-segregated case (j>2 = 0. Rk is the dis- 
tance of the fcth neighbor lattice site from a reference site. 
The dashed lines are 1/v* 



tunneling conductivity. In particular we have evidence 
how this enhancement can be quasi-quantitatively repro- 
duced by using an effective medium theory applied to the 
tunneling resistor network, according to which the effect 
of segregation on the composite microstructure is con- 
tained entirely in the radial distribution function for the 
conducting particles. By using some simple geometrical 
considerations in combination with the effective medium 
approximation we have been able to provide an explicit 
formula for the composite conductivity as a function of 
the conducting filler content and of the degree of segrega- 
tion. When applied to our model of lattice segregation, 
we have demonstrated how the effective medium theory 
closely reproduces the Monte Carlo results for the con- 
ductivity, which can be interpreted as due to a reduction 
of the available lattice sites for placing the conducting 
particles. 
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